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ABSTRACT 


The second-order solution of the problem of the inter- 
action of a train or regular waves with a completely sub- 
merged, horizontal circular cylinder in finite depth water 
is presented for the two-dimensional caSe. The incident 
wave 1S developed as a second-order Stokes! wave by use of 
a perturbation method and the solution of both the first- 
order and second-order scattering potentials is obtained 
numerically using the Green's function approach. The hydro- 
dynamic pressure and resulting first-order and second-order 
SeeeecmcOCtitcients are determined numerically and presented 
for various values of water depth, cylinder depth of 


Submergence, and wave length. 
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P. INTRODUCTION 


The solution to the two-dimensional linear wave/structure 
interaction problem has been well-studied for both the finite 
and infinite depth cases. The fluid motion resulting from 
the interaction of a train of regular waves with a submerged 
homrzZontal circular cylinder in infinite depth water was 
meGEst Stualed by Dean [{Ref. 1]. Ursell [Ref. 7] later studied 
the problem anew placing the solution on a rigorous basis. 
More recently, Ogilvie [Ref. 6] reconsidered the problem. He 
computed the first-order oScillatory forces and second-order 
Steady-state forces for the following cases: (a) the cylinder 
held fixed, (b) the cylinder forced to oscillate sinusoidallv 
iiestil! water, and (c) the neutrally buoyant cylinder allowed 
to respond to wave motion. 

It can easily be shown that the Steady-state part of the 
Second-order forces can be computed from the first-order 
potential alone. Accordingly, Ogilvie did not obtain a 
complete second-order solution; the steady-State forces 
arise from the first-order velocity squared terms in Ber- 
noulli's equation. In addition to the steady-state forces, 
Second-order oscillatory forces also exist which have a 
frequency twice that of the first-order forces and represent 
the subject of the present work. 

This thesis reconsiders again the submerged horizontal 


cylinder problem but the extension is made to include water 
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of finite depth. The cylinder 1s considered to be completely 
submerged and held fixed in a train of regular waves. The 
objective is to determine the ere second-order solution 
and, accordingly, determine the oscillatory second-order 
forces as well as the steady-state second-order forces. 

The problem is treated as a regular perturbation problem 
in the small parameter, incident wave height/cylinder radius. 
Proceeding in this way the incident wave appears as a second- 
order Stoke's wave. The boundary-value problems for both 
the first-order and second-order potentials are established 
and the solutions to both are obtained by use of the Green's 


fumeGtlon method. 


in. 





ieee Orel eA Le DEVELO MENT 


A. FORMULATION OF THE PROBLEM 

A systematic representation of the problem considered is 
mibiecerated inwrig. (1). A rigid right cylinder of radius a, 
Piemercged £o a depth d in water of depth h, is subjected to 
a train of regular waves propagating in the positive x direc- 
t20n. The basic problem is that of calculating the hydro- 
dynamic presSure on the immersed cylinder and the resulting 
force correct to the second-order in wave height of the 
incident wave. 

Pscunmag the tluid £O be arrotational, a velocity 
potential, 6, may be defined as: 


- 


G = Vo(x,y,t) | ay 


where q Genotes the fluid velocity vector. (The barred 
quantities denote dimensional quantities.) Moreover, assuming 
the fluid to be incompressible, the velocity potential must 


Satisfy the Laplace equation, 
— 9) Sees SS = 


woethin the fluid region. 
In addition to the differential equation, Eq. (2), $ 
Must Satisfy certain boundary conditions. In specific, 


these are: 
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Hieszero NOrmMal velocity condition an Poem Ot Lom, 
defined by y = -— h, where h denotes the mean fluid 
depth. 

The zero normal velocity condition on the surface, 
S, (x,Y) —MOemonetie cvlander- 

The kinematic boundary condition on the free surface, 


So COT = in (eee) = 02 


4. The dynamic boundary condition on the free surface. 


These boundary conditions may be stated mathematically as 


follows: 


where n denotes the elevation of the free surface, g denotes 


o—(x,-h,t) = 0 (eS) 
Ve" VS4 = 0 (4) 
. aS, 

° ae — 
q VS. me 0 C5) 


(6) 


the acceleration of gravity, and K denotes the Bernoulli 


eonstant. 


For convenience in carrying out the solution, the 


variables are next made dimensionless using the cylinder 


radius, 


and the wave frequency, o, as follows: 
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x = x/a Y= See ae — da he = ha 
fe H/2a K = K/a v = 0° a/g eS ote (7) 
¢ = ob/ga VS WA 


H denotes the dimensionless wave height, K denotes the 
dimensionless Bernoulli constant, and t denotes the 
dimensionless time. 

Utilizing the parameters defined in Eq. (7), Eqs. (2-6) 
may be rewritten to concisely define the boundary-value 


problem in dimensionless form as: 


rene) = 0 Do ce. f Peie rea Zon (8) 
Xx,~-h,t) = Q 9 

o,, { ) (9} 

o, (%1¥ rt) = 0 on S, (x,y) = 0 (10) 

d(x n,t)ny (x,t) = o, (xen,t) a vn, (x,t) = (11) 


o, (x,n,t) ar slo (x nyt)” ste by (xen t) A] Pen xyt). = h- “Cleey 


B. PERTURBATION EXPANSION 
Expanding $@ and n in terms of a perturbation parameter, 


€, provides a means for converting the nonlinear boundary-value 


ne, 





problem into a series of linear problems. The solution to 
each linear problem may be tractable whereas the nonlinear 
problem is unsolvable. 

Since ¢, n, and K are functions of the small parameter, 


e, they may be written in the power series: 


oe) 


o(x,y,tre) = 2) e 6 (x,y,t) (13) 
n=l : 
nix,tre) = >» eon, (x,t) (14) 
n=l] 
and 
= COR, ee (15) 


ieeas. (11) and (12), » contains n, and, therefore, e 
miemrcitliy. This can be converted to an explicit form in 
e by use of the Taylor Series expansion: 
5m 


S n(x,tre)™ o (xy, t), 
= m a 
m=0 m! oy 


O(x,n,t) = (16) 


0 
The perturbation parameter, ¢€, iS related to the wave height 
in an, as yet, unknown manner. 

Equations (13-16) as well aS appropriate derivatives 
Similar to Eqs. (13-16) may now Ne substituted into the 
boundary-value problem given in Eqs. (8-12) to obtain a 


Series of linear boundary-value problems for oy: bos etc. 
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That is, upon equating coefficients of like powers of €, 
time LOllOowing problems for the first two terms in the 
expansion for $¢ can be precipitated: 


First-order (ce): 


V7, (x,¥,t) = 0 (17) 
by (xh t) = 0 (18) 
dyn ber¥ 1b) = 0 on S, Gx,y) = 3 (19) 
o, y (%10,t) i vn, (%-t) = 19 (20) 
n, (x,t) * o,4 (%,0,t) = 0 (25) 
Second-order lem he 

2 = 
V 5 (x,¥,t) = 0 | (22) 
boy (x,ch,t) = 0 (23) 
bo, (xrY,t) = 0 on S4 Os,y) = 0 (24) 
boy (x,0,t) saz Vo, (x,t) 

(25) 


“Ny (x,t) by Vy x, 0,8) ote O15, (%,0,t) ny, (x,t) 
No (x,t) + $5, (x,0,t) = “Ny (x,t) 7.4 (x,0,¢) - 


an 2 2 | 
Nyy (ert) Oy, (x, 0,t) Fy Poy, (%,0,t) + byy (%,0,t) } + K, 


(26) 


4M 





There iS a striking Similarity between the first-order 
and second-order problems; only the right hand side of the 
free surface boundary conditions differ. In the first-order 
problem the free Surface boundary conditions are homogeneous 
whereas the second-order free surface boundary conditions 
are non-homogeneous functions of the first-order velocity 
potential, first-order free surface elevation, and their 
derivatives. 

The first-order potential function may be represented 
by a function which iS periodic and, therefore, the complex 
pocential, u, (x,y), is defined as: 


7 —-~it 
o, (x,¥,t) = abR,[iu, (x,yle 


] (27) 
Meme R_ denotes the real part, a = 2a ie edenotes tne 


wave length, and b is an unknown real constant. 


oe THE FIRST-ORDER BOUNDARY VALUE PROBLEM 
Since the first-order boundary-value problem is linear, 
the potential 4 may be expressed aS a sum: 


b= op + OF (28) 


I : 
where >) denotes the incident wave potential and $ denotes 
the scattering potential due to the presence of the cylinder. 
In terms of the time independent complex potentials, using 


Eqs. (27-28): 


(bps 





(0 ee a (2.95) 


Remembering that the incident wave potential represents 
only the incoming wave, with all the effects of the cylinder 
represented by the Scattering potential, then the incident 
potential must satisfy the first-order boundary-value problem 
when no rigid body is present. Therefore, $7 satisfies 
Eqs. (17-18) and (20-21), which represents simply the first- 
order boundary-value problem for a train of regular waves. 
The solution to this well-known problem in terms of the 


complex potential is: 


devon. coshfa(yth)] _1ax 
a (30) 


i 
_ ~—S —s : e 
a COsn (ah) 


U 


~ 


The relationship between v and ais defined by 


fee 
y= or = a tanh (ah) er) 
mieten, in dimensional form, using a = ka, and k = 27a/1 is: 
2 a 
o- = gk tanh (kh) (32) 


Since the solution for uy is known, the boundary-value 


Problem for = may now be eStablished. Substituting Eqs. 


i 
Ove 29) ahd (27) into the first-order problem given by 


Bas. (L7-21), and eliminating ni between Eqs. (20) and (21) 
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results in a boundary-value problem for the first-order 


scattering potential as follows: 


Vue (x,y) = 0 (33) 
> (x,-h)_= 0 (34) 
U7 (x, = 


IL lal 


S be, sale . 
Uy, (XY) = sos TERY ["y Sapien) leet in, cosh [a (y+h) |e 


on S, (x,y) = 0 (35) 
Uy (%, 0) - vur(x,0) = 0 (36) 


where ice and ny, denote the x- and y-components, reSpectively, 
Sewcehe unit normal vector, n = in, + jn, directed outward 


imoee the fluid. 


D. SECOND-ORDER BOUNDARY-VALUE PROBLEM 
cite Or the lineariey Of themsecond-Order probien, 
the same procedure as uSed in the first-order problem may 
be applied and leads to the representation of the second-order 


potential as the sum: 
if Ss 
5 = 5 i 5 3 (37) 


where $3 denotes the second-order incident wave potential 


and $5 denotes the Scattering potential. 





Considering the case where there 1S no body present, 
there is no scattering wave and therefore, $7 = $5 = 0 
Additionally, the boundary conditions on the immersed cylinder, 
Egs. (19) and (24) are not applicable. When the substitu- 
fons Of $3 us rT and $5 leh Oo arewmade Intce Bas. (22) 
mmo (25nd (26), along with Eqs. (20), (21), (27) and 
(31), and upon eliminating ny and N> between Eqs. (25) and 
(26), the resulting boundary-value problem for the second- 


order incident wave potential is: 


0° 65 (x, Yt) = 0 (38) 
Ab 
bo. (%.~h,t) = 10 (39) 
boy (%,0,t) + Vo sep (X10, ) ee Sb (a“ a *yR, (gS oe NESE) 
(40) 


The periodic solution to the incident wave boundary-value 


problem specified by Eqs. (38-40) is known and may be written 


as: 


CY) 


$5 = 25 VR, [ius Meyer. (41) 


ia 


where the second-order complex potential, us 


>! 1s given by: 


ae Coed _ Ph. cosh({2a(yth) ] oitax (42) 
2 2a . A 
sinh’ (ah) 
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Bas. (41-42) are familiar in wave theory and exactly the 
same form as that given by Ref. 4 for second-order Stokes' 
waves. 

Having developed a solution for the second-order incident 
Wave potential, bo, the problem for Second-order aGattennns 
potential may be determined. oy = $3 ct oy and 5 = $5 + $5 
are Substituted into Eqs. (22-26), and the previouSly deter- 
mined solutions for the incie€ent wave potentials, Eqs. (30) 
and (42), are utilized. After applying the relationships of 
Bas. (20), (21), (27 )eoe( 29) and) (4im@and, eliminating No 
between Eqs. (25) and (26), the resulting boundary-value 


problem for the Second-order scattering potential becomes: 


Vo5(x,¥,t) = 0 (43) 
Se ee) (44) 
oy xX,7N, 
S 3 ab* Re 
d5,(x,¥,t) aes 7 eh: 7 eee (ny, cosh[2a (yh) 
Sinh (ah) 
oe sinh [2a (yth) }) e42 (@%-t) | (45) 
on S, (x,y) = 0 
S S oa 1s sg noe ae | 
2 2 
St iS S S -i2t 
ieee oi) Oyo 
Us) on y = 0 
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where U(x) 1s a non-periodic function generated by the 
Substitution of oy Teor sow 2zoymancd (260), For brevity, 
this function is not written out since it iS not needed in 
determining the Second-order pressures and forces. 

The boundary-value problem in $5 given in Eqs. (43-46) 
is time dependent according to grees except for the U(x) 
term in Eq. (46). Therefore, the solution for $5 may be 
taken in the form: 

68 = 2 ab*vRjilus(x,y)e 77" + OS (x,y) (47) 
where the last term denotes the time independent portion of 
the complex potential. Separate boundary-value problems for 
us and us Pease mErOmmtiec SUbStrtution of Ho. (4/7) into Eas. 
(43-46). However, aS will be demonsStrated in the pressure 
and force development, only the oot term is required. There- 
fore, the time indepencent part of $5 will be zero, and as such 
will not be further developed; concentration will be directed 
towards the solution for u>. 

When Eq. (47) is substituted into Eqs. (43-46), the 


resulting boundary-value problem for the second-order 


Seattering potential is: 


Vous (x,y) = 0 (48) 
recon = 0 (49) 
Dis eines 


2) 





1 


8 (x,y) = ——b— 
Sinh’ (gh) 


In Samii 2aiyth) | 


+ in, cosh [2a (y+h) } | eee on S, (x,y) =0 (507 


S = S ° = * » a 
Ug, (x,9) 4vu. (x,0) £* (3) (51) 
where 
wale S oO Ls rs IES at 
rs ARGS ant = ca 
Pe GO = ayloMiytiyy * Mitiyy 7 SMrytay * vMiyYayy ©?) 
2 Z 
1s: S S 

- cae ite 7 2S se . Sy 1 y=0 


Again, there is similarity in form between the first-order 
and second-order problems; the only difference in form being 
that the second-order free surface boundary condition, 
Eq. (51), is non-homogeneous. 

peice EC. (51) 1s non-homogeneous, further use of the 


linear superposition theory is made by defining uw 


5 as the 


sums: 


Us = uy ate Us (53) 


O 


* 
where us denotes the homogeneous solution and u> 


Z 
particular solution to the boundary-value problem as stated 


O 


* 
in Eqs. (48-51). More precisely, us and us are defined as 


the solutions to the boundary-value problems obtained from 


the 


the substitution of Eq. (53) into Eqs. (48-51) as follows: 


28 





: * 
RParticuiaxs Solution (us ): 


* 
Vous (x,y) = 0 (54) 
2° a) = 6 55 
Uoy (x,-h) = (55) 
TG Oe neue = en) (56) 
2a 2 d 


Homogeneous Solution: 


Oo 


vous Cains OU oe 
s® 
U5, (xh) = 0 (58) 
So S° 
Woy (x.0) ~ évu. (x70) ="0 (59) 
s° 1 
U5, (XY) nia ne n Simmizal(yrh) | 
San (an) 
. * 
+ in, cosh [2a (yth) } ]e" 7% ws (x9) (60) 


on S, (x,y) = 0 


By the judicious division into homogeneous and particular 


« * e 
solutions, the non-homogeneous problem for us Ceueains 116 


boundary condition on the cylinder surface. This boundary- 


* 
value problem for ub as stated in Eqs. (54-56) is identical 


2 


to that associated with the linear problem for the potential 


Zo 





resulting from a harmonic pressure variation of amplitude 
distribution f*(x) on the free surface in water of depen hk 


O 


The homogeneous boundary-value problem for we , Gefined by 


2 
Eqs. (57-60), is Similar in form to the first-order problem 
given by Eqs. (33-36), the significant difference being the 
term 4yv in place of v which occurs in the free surface boun- 


O 


aia@vecondition. Thus, the method of solution for us will 


be Similar to that of the first-order scattering potential. 


Peer ENE TION OF THE INCIDENT WAVE 

Having developed the first-order and second-order boun- 
dary-value problems, it 1S now appropriate to completely 
define the incident wave height and in the process determine 
expressions for the unknown constants b anda Ko. BOA LAG) 
Beiswee(2i) and (26) for ny and No: respectively, and then 
Substituting the results into the expression for the free 


surface elevation as given by Eq. (14) yields: 
ieee) = — €¢d,, + eae Onna ne Oar At | (61) 
i Je Bite ee ie elas 
il 2 2 ¢) 
- = “- 

For, + bi) + Ky] ues) 
where the potentials are evaluated at y = 0. Eq. (61) is an 
expression of the free surface elevation in terms of the 
femal! potential, and therefore, includes the effects of both 


the incident and the scattering waves. Hence, as it is of 


interest to obtain an expression for the free surface elevation 
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of the incident wave, the scattering potentialS are set to 
zero, 1.e. $ - $s = 0. Evaluating Eq. (61) uSing the known 
solutions for the incident wave potentials, Eqs. (30) and 


(A277, along with Eqs. (27) and (41), then yields: 


; 2 2 Zz. 
et) = ebR, fe? (@Xx7®) y + ‘oe + K, (62) 


‘n ab” cosh (ah) [2 + cosh (2ah) ] R Cotter: 
4 e 


Bah? (ain) 


=f 0(e>) 


where moe) denotes the @Gimensionless Surface elevation 
for the incident wave with no cylinder present. If the x-axis 
1S placed at the mean water line then the constant term must 


vanish and, therefore, 


Z 2 2 
ma (saa) 
Ko = eae rae (63) 


The remaining terms in Eq. (62) are time dependent periodic 
functions, the second-order term at twice the frequency of 
the first-order. 

Defining the dimensionless wave height given in Eq. (7) 
as the difference in elevation between the crest and trough 
of the incident wave, then eb must represent the wave ampli- 


tude. The second-order term in Eq. (62) at twice the funda- 


mental frequency makes equal contributions to surface elevation 





at both the crest and trough, so that it contributes nothing 
to the wave height. Thus, in terms of dimensionless wave 


height, b 1S given by 
bl Sele = AY OBI (64) 


where H denotes the elevation difference between the wave 
Seest and trough, 
ExpresSing the incident wave profile in terms of the 


dimensionless wave height yields: 


pat) = H cos (ax-t) (65) 


Ha cosh (ah) 


cae [2 + cosh(ah) ] cos[2(ax-t) ] 
4 sinh™ (ah) 


=> 
It may be noted that Eq. (65) agrees with the expression 
for the second-order Stokes' wave given by, for example, 


Ref. 4. 


Peewee REOSOURES AND FORCES 
The pressure may be formulated by uSe of Bernoulli's 


equation, arranged as follows: 


2 2 
a a ss “i a _ = 
Bec ©) = po xP LO + oF } egy + pgK (66) 


where P(x,y,t) denotes the dimensional pressure and p denotes 


PMeienrlutad Gensity. By wse of Bas, (7), (13-16), (27), (41), 


SZ 





(47), and (64), Eq. (66) may be reduced to dimensionless forn, 


and carried to the second-order in wave height as follows: 


at 





Peet) — -y — Hak iu,e 1 (67) 
Da 2 
iearal 6V 2 2. ~i2t 
~ Tg [Ret ay Me ea 
2 Ze. 
+ 


v 
ey + buy! + Oy = 2] 


In Eq. (67) p(x,y,t) denotes the dimensionless pressure 


coefficient and is defined as: 


Ps 9 ee) 
pga 


DiGeny pe) a 


The first term in Eq. (67) represents the hydrostatic 
pressure as y is the dimensionless depth beneath the mean 
free surface (y = 0). The second and third terms are harmonic, 
the third term having twice the frequency of the-second and 
representing the harmonic second-order contribution. The 
remaining terms in Eq. (67) are independent of time and 
provide the time-average or steady-state force components. 

Expressing the components of the wave force vectors in 
terms of integrals of the pressure over the cylinder surface 
area, the dimensionless force coefficients may be written 


ass 


C, (t) = -_fetxry ting ac, i De Use) 
i 





where the dimensionless force coefficients in the x and y 


directions are defined, respectively, as: 








F(t) 

Ce) at (70) 
pga 
Fe) 

Cy (t) = i; (7a) 
pga 


with Eee) and ae denoting the force components. Addi- 
mironally, aC, denotes a dimensionless differential arc 
length along the circumference of the cylinder. 


Applying Eq. (67) to Eq. (69) yields: 


eee _ = ie e 
C; (t) = y frye, Ha fRetuye Jn, ac, (72) 
c c 
aR 
2 2 : 
2,;a 6V 2 De = 1 2 
oT Mires _ SP a2 ex Wye 
] 
2 2 2 3 
~ + 
+ Jus! a1, + 2 1]n,dc, + 0(H™) 


As the first term in Eq. (72) results from the hydrostatic 
pressure on the cylinder, it represents simply the buoyancy 
force. Omitting the hydrostatic force, the hydrodynamic 


force may be written as: 
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7 _ 2 
Cc, (t) = HF cos (64; t) + H (73) 


Ie. 


2 7 Ss ane 
+ H [Fo 5 cos(62i-2t) + Foy! ae (Use) 


where the first-order, second-order, and steady-state force 
coefficients and phase shift angles are defined by comparison 


of Eqs. (72) and (73) as follows: 


Be § 
Fie =a un, dc, (74) 
C 
1 
16 2 
7A ss ae! f& . Zo 2 
F.;¢ = av (——u, My) ie or (75) 
Cc 
1 
2 2 2 2 
SS a f Vv 
— = —— Groweier fae | see nae (76) 
Z1 Ay o 1x ly a eer 


The dimensionless force coefficients, Pia and Poy are real. 
The first term in Eq. (73) represents the first-order 
forces of the fundamental frequency, o. The periodic portion 

Gf the second term in Eq. (73) represents the second-order 

contribution to the force at twice the fundamental frequency. 

The last second-order term represents the steady-state or 

time-independent contribution, sometimes called drift force. 
Evaluation of the force coefficients defined by Eqs. | 

(74-76) is the main objective of this thesis. Therefore, 

it 1S necessary to evaluate the first-order and second-order 


complex potentials, Uy and u as well as derivatives of Uy 


go 





on the surface of the cylinder. Additionally, evaluation 


of u, and its derivatives on the mean free surface (y = 0) 


18 


is required. 
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fii; METHOD OF SOLUTION AND NUMERICAL PROCEDURES 


Pee oOLUTION OF THE FIRST-ORDER PROBLEM 

The use of a Green's function to express the solution 
to the first-order boundary-value problem was first formu- 
lated by John [Ref. 5], and applied to submerged ellipsoids 
by Garrison and Rao [Ref. 3]. This method is considered to 
be applicable, in principle, to arbitrary shapes and is 
mathematically the most Straightforward. The Green's func- 
tion, G, of unit strength which satisfies Eq. (33) as well 


eomeeie boundary conditions, Eqs. (34) and (36) is given by: 


a 7 cosh[u(yth) |coshtiy (nth) ] 
GCx,yr iene) Pia 2 ae eae apy _f [eos tutynd leper usinh uh) 
7 oe HA ik. oe Ale = 
Satan uh cos|x-&| dy | (77) 
: At 
- § ZaNFSINATZapRT Cosh La, (Y#H) 1eosh fay (n#h) Joos fay |e-6 | 
where: 
L 
c= [(x-€)7 + (y-n) 71° (78) 
i 
r= [(x-£)*% + (ytn) 71? (79) 
The symbol, ay 1s defined in terms of h and v as the 


solution to the equation: 


On 






F(a,,h,v) = 0 (80) 


where F is defined by: 


F(a, ,h,v) =~ ay tanh (a,h) - Vv (81) 


Comparison of Eqs. (31), (80) and (81) demonstrates that ay 
1s clearly the equivalent of a. However, this will not be 
the case for the second-order problem, requiring the use of 
separate notation. In Eq. (77), PV denotes the principal 
value of the integral. 


This Green's function was also given in Series form by 


JOitmtRein 5} as: 


Zone? 
00 (uty )cosfu, (yth) ] =| oaks 
Slee yt ee v) = 27 2 Se cosy, (nth) Je 
k=1 Vu Thuy ( Wty ) 
ia,|x-€| 
4ncosh [a, (yth) ]cosh[a, (nth) Je 
Sr Vor) 
2a,h 3 sinh (2a,h) 


where ay 1s as defined in Eq. (80) and the quantities Hy 


are defined as the real positive roots of the equation: 


K (uy, hv) = 0 (83) 


mhnere the function K is further defined as: 
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RCH) eee one ( yee (84) 


Following the Green's function method of solution, us 
femwratten as the integral over the cylinder are length, 


Cc as: 


cag 


ve (X,Y) = = fejemocevieniv) dc, (85) 


al 


where (&,n) denotes points on the immersed Surface, £,(&,n) 
denotes the unknown source strength, and dc, = dc, /a denotes 
the differential arc length on the cylinder surface, mace 
dimensionless with the characteristic length, a. 

Although the solution to the first-order boundary-value 
problem as stated in Eqs. (33-36) is given by Eq. (85), the 
source strength function, £, (Eon), must be determined in 
order to evaluate the potential. From potential theory, 
the normal derivative of the potential, us (X,Y) , on the 
surface of the cylinder is: 


S ai 1 ae 
Us, (*-¥) ae) £, (x,y) te oe _ Jiri, GoysEnivracy (86) 


al 
where Gas the normal derivative of G, may be determined by 
differentiation of either Eq. (77) or (82) in a straight- 
forward manner. 
Applying the final boundary condition, Eq. (35), leads 
to an integral equation which may be solved for f 


lit 


BY 





a) th Be fey tems, tvtime) dc, (87) 


1 


Nike 


- iL : : 1ax 
= soci tany |PySinhla(y+h) | Ta in cosh [a(y+h) } Je 


The solution for fy from Eq. (87) may then be used in Eq. 


(85) to determine the potential, Hee 


Bee soOLUTION OF THE SECOND-ORDER PROBLEM 

The Similarity in form of the boundary-value problem 
for the homogeneous part of the Second-order potential, 
Boas.) (57-60), to the first-order potential is now utilized. 


Since the only significant difference Occurs in Eq. (59), 
O 


where v 1S replaced by 4v, we may represent us in a form 
Similar to Eq. (85) as: 

eve 1 

uy (x,y) = oe _Jigteim tase vite) ac, (88) 


al 


where G(x,y;&,n;4v) is defined by replacing v by 4v in Eqs. 


(77) and (82). Therefore, a. is defined by: 


2 


F(a,,h,4v) = 0 (89) 


and a. iS replaced by a. in Eqs. (77) and (82) also. In the 


il 
Seese Of Eq. (82), 


2 


Hy is defined as the real positive roots of: 
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K(u,,h,4v) = 0 (90) 


The source strength function, £i(E,n), appearing in 
Eq. (88) is again determined by applying the kinematic 
boundary condition on the cylinder surface as given by Eq. 
(59). The integral equation for £. may then be given by: 


= f(x,y) + Ae fegteoms, GeavsEonzav dc, ene) 


il 


— + in, cosh [2a (yth) ] 


i2ax S* 
Sinh (ah) 


on S, (x,y) = Q 


os 
However, tO SOlve Eq. (91) for fos Ei eotecrsOlLuLIOn Lon U5 
* 
must be obtained and tee evaluated on the cylinder surface. 
, 7 
tiie, at this point the solution to Se is sought. 


2 


x : ‘ 
The complex potential ce is defined as the solution to 


2 
the boundary-value problem, Eqs. (54-56). The form of this 
problem iS recognized as being equivalent to that of fluid 
motion produced by a periodic pressure distribution (as 
specified by f£*(&)) with frequency 20 on the free surface. No 
cylinder is considered to exist in the fluid region. 

The solution to this problem may be formulated as an 


integral over the free surface extending from negative 


mofinity tO positive infinity. Using the present 
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* * * S* 
dimensionless representation, U5 becomes: 


3.8 ge 7 ; ; 
Us ey) a OT jf: Ue ue (x YY eae 4) dé (92) 
, Cs 


The Green's function, G*, for the problem is given in Ref. 8. 


Upon transforming the solution to the non-dimensional form: 


~ cosh{u(yth) Jcos[u(x-£&) Jau 
—0o 4vcosh(uh) - psinh (ph) 


Ger, y;£,0;4v) = =-2PV 


4mcosh [ah] cosh [a, (yth) ]cos[a, (x-&) ] 
een (9.35) 


2a 5h -|- sinh (2a.h) 


where a. is defined by: 


2 


F(a,,h,4v) = 0 (94) 


o 


* 


The normal derivative of us is required in order to 


completely define the kinematic boundary CcOncHe1on om the 


cylinder as given in Eq. (60). From Eq. (92), this becomes: 
ue (x,y) = te ER (EG. (xyy7E,0;4v) aE (95) 
2n ‘*s 27 Ay eee 
2 


p@emaerivative of G* in the direction normal to the cylinder 
surface may be obtained by differentiation of Eq. (93) ina 
Straightforward manner. Thus uSing the values for £f*(&), 


as defined in terms of the first-order solution only, as 
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. Ss So: 
given in Eqs. (52) and (56), u, and uz, 
O 
* 


’ US may then be determined to complete the 


may be determined. 


; S 
(oissal ale U5, 


second-order boundary-value problem solution. 


C. NUMERICAL METHODS 
Determination of the forces on the cylinder surface 


requires solving for the potentials u us, and the derivatives 


ee 
Ot Uy at points on the surface, as shown in Eq. (72). Thus, 
the problem is now one of solving for both the first-order 


and the second-order scattering potentiais, which requires 


solving for the source strength functions f, and f as well 


iL 2s 


as the free surface pressure distribution source strength 


function, £*. In view of the complexity of the equations 


Coad 


. : ae . at er eae 
Mmampesco mature! £O attempt <= numerical sclution. 


The first~order scattering potential, ae aS given by 
Eq. (85) is dependent upon the first-order source strength 


momecton, f£ Thus, it iS necessary to solve Eq. (87) for 


1 
fy to determine ae A numerical solution may be developed 
by dividing the cylinder surface into elements of length 

A6é = 271/m, with the center of each element assigned an index. 
Since £, (én) 1s recognized eas a well-behaved function for a 


Smooth surface cylinder, it iS appropriate to define the 


following: 


=) 


. Jog taxpegreoni dc, (96) 
a) 


O54 (v) =. 
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sua | 
h; = cosh(ah) pry Oy r¥,) sinhfa(y,+h) ] 7) 
Lax, 
+ in, (x, /y,) cosh[a(y,+h)] | e 


and thus Eq. (87) may be approximated by the complex matrix 


equation: 


£43 -- O44 (vf). = 2h) - Deis aa al be a (99) 
Upon evaluation of O34} uSing Eq. (96), the inversion of 
pepe 99) may be carried out on the digital computer to 
determine £5 at each nodal point on the surface of the 
cylinder. 

For the purpose of evaluating a and 8B, either Eq. (77) 
or Eq. (82) may be used. For evaluation of a and B when 1=j) 
Eq. (77) must be usSed to allow separate treatment of the 
Moeerithmic Singularity as r approaches 0. The difficulty 
encountered with the logarithmic Singularity may be overcome 
by carrying out its integration analytically. Garrison 
[Ref. 2] showed for the calculation O as Giacene Comer lou 
Elen Or the normal derivative of the ln (r) in Eq. (77) 
integrated over the singular element of A@ is A6/2, not only 


for the nodal point (x. /Y,) — (er ioeeoue tom all nodal points 
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(x5 7Y,) on’ the cylinder surface. Additionally, to calculate 
Bs Garrison developed a numerical approximation for the 
integral of the in (r) over the Singular element ene fA\lele 

A second Singularity arises in oe evaluation of the 
ieimece integral Occurring in Eq. (77). The first term in 
the integrand 1s Singular at up = a. Since the numerator 
approaches zero like (u-a) near a, Garrison [Ref. 2] Sub- 
mateeed 1/({\i-a) from the Singular term in the range 0 < u < 2a 
and added the -contribution of the integral of 1/(\-a) which 
in this case waS zero. This converted the Singular integrand 
to a regular function. Applying this technique to numeri- 
cally evaluate the resulting integral between 0 and 2a, 
Simpson's three-eights rule is used with an odd number of 
Subdivisions, thus ensuring that the point pw = a is not 
encountered. 


With f, determined from the inversion of Eq. (99), the 


i 
first-order potential and its derivatives may be evaluated 
on the cylinder surface. Replacing the Surface integrals 


With summations, Eq. (85) and its derivatives may be written 


ass 


oe = 8.5 (v) fy, al = 1 Die Sipe cael (100) 

i = pee. (101) 
iat alg sls} 

noe Ge ()) 58 (102) 
ly Vas) 13 
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ae ete. denote functions evaluated at™the 


S 
where u Lyi’ 


i aglag 
ae nodal point on the cylinder Surface. The complex 


Matrices in Eqs. (100-102) are defined as follows: 


ee : ; Sanhe aes 
B.4() moe ovr, vgreiniviac, US a ctl ro 2) ears (2 ies Gl 8)S)) 
KC 
zy 
B (v) = = GG, ena) ae (104) 
sol) 27 n > dealt eas rane 1 
1G 


ee i ; 
a? Shor Js, x. 7Ys5S,niv) dC, (105) 
AC, 5 


The first-order incident potential and its derivatives may 
be evaluated directly at each nodal point after differentia- 
tion of Eq. (30) as needed, and combined with the scattering 
potential and its respective Gerivatives. 

To solve the second-order problem, the function f¥*(&) 
required in Eq. (95) and as defined by Eq. (52) must be 
evaluated numerically. Dividing the mean free surface 
(y = 0) in the vicinity of the cylinder into n equal incre- 
ments, £* 1s evaluated at each nodal point. The numerical 
solution uses again Eqs. (100-105) with eq Omeote at he 
purpose of evaluating ee and its derivatives at the mean 
free surface nodal points. The Migidenierirc t-Oncer Potential 
and its derivatives are evaluated directly at each nodal 
point, and combined with the scattering potential and its 


derivatives to solve Eq. (52) for f*., 
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Having determined f* at 
surface, the boundary-value 
(54-56) may be solved, l1.e. 
on the cylinder by solution 
tavely. 


complex matrix form yields: 


aries: nodal points on the free 


* 
problem for ue given by Eqs. 


2 
ot Sie l 
Us and Voy may be evaluated 
Of Eqs. (92) and (95) respec- 


Expressing the integrals as numerical Summations in 


= : 
u, . = £*a%. (4v) 1 =e 2, Pall (106) 
Soya a ee a += 1,2, i 
S* _ £xB*, (Ay) (107) 
an 5 uy 
where 
4 (4v) = = Gr ( ; Oy) ela (108) 
iy Vv = Xi 7Ysi6sNi V 
NG 
2) 
1 
* a oan * s e ; 
BF, (4v) aT Jor cxgyg:tensav ells (109) 
INO 
2) 
S* Ss ; 
imosee using £*, both u ana u may be directly evaluated 


2 


2n 


at each cylinder surface nodal point. 


To sOlve the boundary-value problem for the second part 


of the second-order scattering potential, us 


by Eqs. (57-60), Eq. 


O 


(88) for 


evaluate Eq. , 


(88) must be evaluated. 


O 
, aS Specified 


However, to 


mete soe cond.Ordemm@source Strength 


Punction, f., must be evaluated by numerically solving 
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the integral equation, Eq. (91). Replacing the integral 


equation with a complex matrix Summation: 


F.. + £ 


oe A ae 12 (110) 


where a; 4 (4v) is @efinede by EG. (96), 


to = Eo (ey Ys) (aed 1) 


k. = Peete ae) [n 


Rez. ) Santi | 2aty eenye) (aneZ:) 
z Pr es ies Zi 


¥ 


i2ax. Gx 
+ in, (x, ,¥,)cosh[2a(y, +h) ie us, (X14) 


Once Sig BERD, is evaluated, Eq. (110) may be inverted to obtain 


the second-order source strength function, f., at each nodal 


powmmearon the cylinder surface. Stating Eq. (88) in summation 
FOLmM: 
5° 
— a ; 5 —_ 
OH £55815 (4? 1,J U2 op ee (713) 


where Bs 4 (4v) TSmdcminecdeuayere. (hos). Usang Ba, (113), 
O 
S 


u, May be determined at each nodal point on the cylinder 
* 
Surface, then combined with us and the second-order inci- 
. I 
dent potential, Us, to evaluate the total second-order 


potential, Uo- 
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Since the first-order potential, Use and its derivatives, 


Uy and Say) 


now known at each nodal point on the surface, the first- 


, aS well as the second-order potential, U5, are 


order, second-order and steady-state force coefficients and 
phase shift angles may be determined using Eqs. (74-76). As 
shown earlier, the integrals may be replaced by Summations, 
uSing nodal point values and the arc length increment, Aé. 
Once each integral is evaluated, then the force coefficient 
becomes the absolute value Or modulus of the complex integral 
result and the phase shift angle is the angle whoSe tangent 
is the imaginary part over the real part of the complex 


integral result. 


Pee COMPUTER SOLUTION 

In order to utilize the method of solution described 
in Sections B and C, a computer program was developed to 
Carry Out the indicated calculations. Although accuracy was 
a prime consideration, an equally important requirement was 
to minimize computer run time. To achieve this, every attempt 
was made to utilize symmetry and to generate certain constants ~ 
and matrices to eliminate redundant computations. 

The program consists of a main program that flows in the 
Order of computation presented in Sections B and C, along 
with four subroutines. Two subroutines to solve the first- 
Order Green's function, uSing both forms of the function, 
Eqs. (77) and (82), were developed by Garrison [Ref. 2]. 


These subroutines, GREEN and GREENS respectively, have been 


49 





modified to include calculation of the first-order Scattering 
potential derivatives, evaluation of the first-order scattering 
potential and its derivatives on the mean free Surface, y = 0, 
and evaluation of the modified Green's function, G*, for 


* 


* 
pieeaetermination of ue and ne 


; on at each cylinder nodal 


Peint. 

For elements of the a and 8 matrices corresponding to 
small values of the parameter (x-&), GREEN is used while for 
larger values of (x-&), GREENS is used. With the exception 
of the diagonal elements in Eqs. (96) and (103-105) and a few 
surface nodal points in Eqs. (108-109), the majority of 
elements of a and 8 are calculated by GREENS. This is most 
fortunate as the series form converges rapidly, requiring much 
less computer time than the integral form. 

Subroutine GEODAT reads the input geometrical data, 
generates the matrices hs and k. as defined by Eqs. (97) and 
(112) respectively, and calculates certain geometrical 
parameters and matrices for repeated uSe in GREEN and GREENS. 
Subroutine COMAT inverts the complex matrix equations and 
thus is used to invert Eqs. (99) and (110) to determine 
fy and f. respectively. 

A cross-reference between text and computer program 


nomenclature is given in Table I. 
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PAD: 1: 


text 


Computer Program — Text Symbol Cross~Reference 


Computer 
Program 


A 
D 
) GEL) 
ae, L) 
ase Ir), 


Cl) 
Gila) 
C2q) 
C22) 
€3(1) 
G52) 
Cig, GLI 
GId ,GIJEXT 


GNIJ,GNJI 
GXIJ,GXJI 
GYIJ,GYJI 
GYY 
H 
HH (I,1) 
PKG, 1s) 
NPTS 


NSPTS 


Text 


J 


1 (surface) 


u 
ae (surface) 
bx 


al 
Ny (surface) 


we 
lyy 


(Surface) 


us (surface) 


. 
Wy (surface) 
5 (surface) 

ly 

Ae (surface) 
lyy 

bay, 
D4 
y 
a 
B 
8. 
By 


NOMCGY lamider ) 


AE (surface) 


Computer 
Program 


ANX (I) 
ANY (T) 
Ui) 

Bix (hy) 
ey Ga) 


UiLrs 


iso 


LOLS SONG 


eZee) 
ser) 
Y¥ (I) 
ALPHA (I,J): 
BETA (I,J) 
BETAX (I,J) 
BETAY (I,J) 
DELTHE 


DELX 






Cempuver CGMmpuUEe rs 


eae Program oo Program 
644 PHASE1 (1) u -  _AMU 
S45 PHASE] (2) Uy (v) AMU (K) 
Soy PHASE2 (1) Hy (4v) AMU4 (K) 
S55 PHASE2 (2) Vv ANU 
E X (J) 4y ANU4 
n ¥ (J) T Pea 





IV. DESCUSSION AND RESULTS 


A. SELECTION OF COMPUTER PROGRAM PARAMETERS 

The computer program was developed on the IBM-360 
computer uSing FORTRAN H. All Subroutines were compiled on 
DATA CELL to minimize compile time and core size. In order 
to maximize accuracy while minimizing computational time, 
two key computer parameters were evaluated to determine 
optimum values; cylinder nodal points and free surface model 
points. 

Using the first-order solution only, the number of nodal 
points on the cylinder surface was varied from 12 to 60 and 
compared with the results determined by Garrison [Ref. 2]. 
The minimum number of cylinder surface nodal points that 
provided accuracy within one percent of: those Obtained by 
Garrison uSing 60 points, was 24. Accordingly, the good 
correlation from 60 points down to 24 points led to the 
selection of 24 nodal points on the cylinder surface. 

Since the free surface integral was to be carried out 
from -- to tm, the next step was to establish the finite 
interval of convergence On the surface as well as the sub- 
division size. After the outer limits of the free surface 
integral were determined, the subdivision size was eee 
from 4 to 128 subdivisions per second-order wave length 


holding the total interval constant. Above 64 subdivisions 
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Pee eecsules did not vary more than One percent, leading to 
the selection of 64 subdivisions per second-order wave length 
On the free surface. 

Additionally, the total interval, 2, waS varied from 
one to eight wave lengths in both the positive and negative x 
directions. The second-order force coefficients and their 
respective phase shift angles converged to a small value 
varying periodically. Convergence to this limit cycle 
Occurred between the first and second wave length out from 
the center of the cylinder. Therefore, a total interval of 
from -2 to +2 second-order wave lengths was chosen as adequate 
for generation of the numerical results. 

Although the values of the Second-order force coefficients 
and their respective phase shift angles tended to converge 
with increasing the limits of the free surface integration, 
there remained, in general, the small limit cycle as noted 
above. The amplitude of the cycle waS a Eene cion of the 
parameter, a, and was of Significant magnitude only for 
values of a less than 0.25. The effects of a on the limit 
cycle for one of the second-order parameters, horizontal force 
coefficient, is demonstrated in Fig. (2). This represents a 
plot of percent variation of the force coefficient from the 
mean versus the Surface interval size, A/(L/2), corresponding 
to values of a of 0.25, Ue Omand 075. Figure (2) 1s 
representative of the behavior of all second-order force 


coefficients and phase shift angles and demonstrates the 
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increase in the amplitude of the variation with a decrease 
i a. The amplitude of the periodic variation ranges from 
ite, tO Flfty percent at a wave number OD Ol omanea cepenl or 
d= 1.4, to less than two percent for all values of wave 


number above 0.25. 


Bs RANGE OF APPLICABILITY 

In a complex computer program of the type developed to 
carry out the present numerical scheme, certain limits with 
respect to numerical stability, etc. show up. These computing 
limits arise for various reasons such as overflows caused by 
computing hyperbolic functions of very large numbers, 
numerical convergence instabilities, etc. While no attempt 
was made to investigate each of these numerical problems, a 
list of the iMiteasiens of the particular program developed 


in this thesis to carry out the computations are as follows: 


minimum a — dependent upon acceptable error, but 0.25 
or less 


maximum a — deep water condition reached 
minimum h — 3 


maximum h — deep water condition reached for a > 0.25, 
au 2 0 GhOimea 8 0.25 


Minimum d — from 1.4 at h equal to or less than 5 down 
GOu2 .Seae Do = 20 


minimum SMIN — 0.12 when cylinder depth is other than 
near the free surface 


maximum SMIN — 0.30 when cylinder depth is near the 
free surface 
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ee RESULTS 

A representative set of results for the first-order and 
second-order force coefficients have been generated for a 
water depth of h = 5.0, Submergence depths of d = 1.5, 2.0, 
meee and 3.0, with a ranging from 0.15 to 1.2. All of the 
numerical results presented herein were based on 24 sub- 
divisions on the circular cylinder and 64 subdivisions per 
second-order wave length on the free surface integration. 
The surface integral was carried out from x = -L to +L Since 
this was found to be adequate for convergence. Since the 
second-order wave length is half the first-order wave length, 
L, the total interval, is equal to four second-order wave 


lengths. 
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The first-order hori 
are presented in Figs. (3) and (4), respectively. It may be 
noted that, in general, the forces decrease with depth of 
submergence according to expectation since the wave action 
dies out with depth. 

The second-order horizontal and vertical force coefficients 
are presented in Figs. (5) and (6), respectively. Generally, 
the results show the second-order effect to be relatively more 
miteetant as the cylinder approaches the free surface. It 
might be suspected from this that the second-order contribution 
would be much more Significant in the case of surface piercing 


Or ftioating bodies. 
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The ero and vertical steady-state force coefficients 
are oe in Figs. (7) and (8), respectively. These results 
show that the horizontal steady-state force coefficient is very 
small in general. The horizontal force can be shown to be 
proportional to the momentum flux of the reflected wave and 
the reflected wave is small. In fact, in the case of infinite 
water depth, Dean [Ref. 1] showed that the reflection was 
exactly zero and accordingly the steady-state horizontal force 
1s zero. 

It may be noted that the steady-state vertical force is 
positive. This results from the fact that the velocities are 
larger on the top of the cylinder and accordingly, the 
pressures reduced. 

The pnase shift angles of the first-order and second-order 
forces are shown in Figs. (9) - (12). 

To demonstrate the effect that the second-order terms have 
on the forces on the cylinder and their respective phase shift 
angles, a comparison was made between the first- and second- 
order results. Figures (13) - (18) are plots of the horizontal 
and vertical dimensionless forces on the cylinder versus time 
Over one complete wave cycle for both the first-order solution 
alone and the combined first-order and second-order effects. 
Additionally, a plot of the incident wave is included to 
demonstrate the phase shift involved (the amplitude of the 
incident wave has no significance). A mean water depth, h, of 
5.0, a cylinder depth, d, of 1.5 and a wave height, H, of 
0.5 were used, with the wave number, a, asSigned three values, 


tezo, 0.5 and 1.0. 
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eC ORC LON S 


A computer program has been developed to carry out the 
numerical solution of the second-order wave — cylinder 
interaction problem. The first-order, second-order, and 
steady-state force coefficients were determined for the 
submerged horizontal cylinder. 

An approximate method Poye~ ela ne with the second-order, 
nonhomogeneous free surface boundary condition was developed 
which appears to converge except at small values of a, il.e., 


large wave lengths. 
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